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Bell type inequalities are used to test local realism against quantum theory. In this paper, we 
consider a two party system with two settings and two possible outcomes on each side, and derive 
equaUties in local theories which are violated by quantum theory by a factor of 1.522 tolerating 
0.586 fraction of white noise admixture which is twice that of the previous results. 
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I. INTRODUCTION 

The idea of local realism opened its way into quantum theory by the work of Einstein, Podolsky and Rosen which 
is well known as EPR paradox [1], though it is not really a paradox, because it has been argued recently that the 
assumptions of EPR were wrong [2]. However it was left to John S. Bell who derived an inequality based on local 
theories and proved that it was violated by statistical predictions of quantum theory [3, 4]. Since then many attempts 
have been made to derive Bell type inequalities which are violated by a stronger factor, so that it can be tested in 
the real experiments in which errors are inevitable. Among these are Clauser, Horne, Shimony and Holt inequality, 
the so called CHSH inequality [5, 6]. 

Nowadays, with the growing power of computers, the numerical methods have attracted attentions for constructing 
these inequalities as much as possible [7-9] though analytical approaches are in progress too [10]; with the hope that 
some of which could be violated by a stronger factor. 

In this paper we first introduce a simple but general method for constructing Bell expressions which can be applied 
to two party experiments each measuring two observables with different outputs. Generalization to more parties 
and/or more observables is straightforward. Then in a special case of a two party system with two settings and 
two possible outcomes on each side, we use these expressions and show that according to local theories there exist 
equalities which are violated by quantum theory by a stronger factor than Bell type inequalities. Meanwhile we derive 
an inequality in this case which is violated by a factor of 1.621 tolerating 0.293 fraction of white noise admixture. 
Although in this case the tolerance is the same as those derived previously in the literature the range of violation is 
greater. 

II. DERIVATION OF BELL TYPE EXPRESSIONS 

Let's consider a two party system with no known interaction between its parts. Suppose the left party, say A 
performs two possible measurements a with outcomes i € {0, • • • , m — 1} and a' with outcomes i' G {0, • • • , to' — 1}. 
Similarly the right party, say B, performs two possible measurements b with outcomes j € {0, • • • , n — 1} and b' 
with outcomes j' € {0, • • • ,n' — 1}. For simplicity, from now on we use unprimed/primed variables and indexes 
for the first /second measurement for each party whenever applicable and label such a system as mn(i)m'n' . As a 

consequence of locality, for measurements which are not simultaneous, there exists a probability q"Jli,i defined as the 
probability that measurements of a results i, a' results i' , b results j and b' results j'. The total number of q's, Nq, is 
Nq = mxm' xnxn' . As g's are all exclusive and cover all probable events, with the assumption that the probability 
distribution for each measurement is normalized we have 

E €jC = 1- (1) 

These q's are the building blocks of all possible Bell expressions for local theories. 
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If P^l denotes the probability that in a particular experiment, A measures a with outcome i and B measures b with 
outcome j, we can write 

Pai = J2CJw, and SO on... (2) 

i'J' 

There are totally Np number of P's, where Np = (m + m'){n + n'). Eq. (2) in matrix form would become 

P = MQ (3) 

where P is an Np x 1 column matrix. Q is an Nq x 1 column matrix and M is the conversion matrix with dimension 
NpxNg. These P's are not all independent. The rank of the matrix M, denoted by Nj, is the number of independent 
P's. For 22022, 23022 and 23023 the rank of the conversion matrix, M, is 9, 12 and 16 respectively. Please note 
that these do not agree with the results obtained in [11] in which the number of independent P's for 22022 case is 
predicted to be 8, and for 23022 case it is 11 (for m = 2, m' = 3, n = 2 and n' = 2) or 14 (for m = 2, m' = 2, n = 2 
and n' = 3) and for 23 23 case it is 20 (for m = 2, m' = 3, n = 2 and n' = 3) or 19 (for m = 3, m' = 2, n = 2 and 
n' = 3). 

To prove these numerical results analytically we search for constraints which justify the rank of the matrix M. Two 
groups of constrains are directly derived from eqs. (1) and (2). 

The first is normalization of probability distribution in a measurement. That is 

E^a^ = 1' (4) 

E^i' = 1' (5) 

T.Pa'i = 1, (6) 

i'J 

E^'i' = 1- (7) 

However, with the assumption of CH (see [6]) we can write 

pij pi' J' _ pi' 3 pij' _ Ji' n' (o\ 
^ab^a'V — ^a'b^ab' ~ %a'bb'^ \°) 

and from eq. (1) we have 

E Pab^u'V — E laJw — (9) 
i,i',3,j' i,i',j,j' 



E Pa'ip^' = E = 1- (10) 

i,i'J,j' i,i',jj' 

So, for instance, eqs. (4), (5) and (9) can be considered as independent and these equations impose only 3 independent 
constraints on P's. Please note that in this case eq. (10) is automatically satisfied and eq. (6) can be derived from 
eqs. (5) and (10). Similarly eq. (7) can be derived from eqs. (4) and (9). 
The second group of constraints derived from eqs. (1) and (2) are 

(11) 
(12) 

(13) 
(14) 



y^.pgb 

j 


i' 


/ j^a'b 
3 


— / j^a'b' 
3' 


E Pab 

i 


~ / . ^a'bi 
i' 




~ / .^a'b' 



3 



which imply no-signaling. The total minibcr of constraints from the above equations add up to m + m' + n + n' . 
However in each group of the above constraints, one of them is a linear combination of the other, for example from 
eqs. (4) and (5) one can write 



ij' 



and from eq. (11) we have 



Subtracting eq. (16) from eq. (15) results in 



E 



(15) 



(16) 



which is one of the constraint in eq. 
(m - 1) + (m' - 1) + (n - 1) + (n' - 1). 
The total number of constraints, Nc, given by eqs. (4 - 14) would become 



E^ai = E^i'- (17) 

i i' 

(11). So the constraints from no-signaling would be 



Nc = m + m' + n + n' 
So there are only Nj number of independent P's, where 



1. 



Ni = Np- Nc 

= (m + m'){n + n') — {m + m' + n + n' — 1). 



(18) 



(19) 



which agrees with the numerical results mentioned before and is symmetric with respect to interchanging m and m' 
(and of course n and n'). Clearly without CH assumption i.e. eq. (8) we couldn't get this result. We would like to 
emphasize that this assumption is only used here to prove eq. (19). It has nothing to do with the rest of this paper, 
especially the main results which will be obtained later. 

Generally if B is a Bell expression for local theories and —d < B < c with non-negative c and d, then B must satisfy 



— E -^stklPst 

s,t.k.l 

= E (/^'^ 



\ a jj 

— l^ii'jj'jQaa'bb" 



IJ-i: 



W jj' 



(20) 



where c{d) is the greatest of non- negative real numbers /x's(j/'s). 

We have solved eq. (3) numerically for 22(8)22, 23022 and 23023 cases to find all possible expressions that satisfy 

eq. (20) and the complement of each one, i.e. all pair of expressions whose sum add up to 1. However we do not 
discuss it here because the results that we are going to use in the next section can be tested directly and easily. 

III. VIOLATION OF EQUALITIES AND INEQUALITIES 

Using the numerical method mentioned in the previous section we have found two expressions in 22022 case which 
are complement of each other. These are 



I plO _ pll I pll I pOl I _ „0001 , 0011 , Olio , 0111 , 1000 , 1001 , 1100 , 1110 
Kll -f^l2 + -f^21 T -f^22 I — 91212 + 91212 + 91212 + 9l212 + 91212 + 9l212 + 9l212 + 9l212) 



(21) 



and 



pUl I pOO I pOi I p 
-f^l 1 "T -'^12 "T -'21 T -f^; 



„oooo , „ooio , „oioo , „oioi , „1010 , „1011 I „1101 , „1111 
9l212 "T 9l212 "T 91212 9l212 "r 9l212 9l212 "r 9l212 9l212 • 



(22) 



Please note that all g's are non-negative numbers. Adding these two equations and then using eq. (1), would give 

Pll + Pll\ + I - ^^I°i' + ^1°2° + ^2°i' + ^2^2^ = 1- (23) 



\Pll 



pll 
-'12 
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Sec Appendix A for a direct proof of eqs. (21), (22) and the above equality. (Though the equality (23) can also be 
derived analytically from eqs. (4 - 14) easily.) Of course, due to symmetry, there are other equations of this type as 
well which we do not mention here for brevity. 

Now we use the same experiment used by CH [6] to show that equality (23) is violated by quantum theory. Consider 
a two photon system in the state 




(24) 



where one photon moves to the left in the +z direction and the other moves to the right in — z direction. The 
projection of I'i'n) on directions u(^) = cos^i + sin^j + Ok (for the left photon) and v((/)) = cos(/)i + sin(/)j + Ok (for 
the right photon) is: 



cos^{0) cos(6') sin(6') 

\^{e,4>)) = { cos{e)sm{e) sm'^{e) o 





cos^((/)) cos(^)sin(^) 
cos(^) sin((/)) sm^{4>) 




l*o). 



So 



|*(^,.^)) = -^cos(^- 




(25) 



If for example the left photon is not detected in the u(^) direction then it must be detected in the direction u(0+7r/2) = 
— sin^i + cos^j + Ok which is perpendicular to u(^). So defining 7 = ^ + 7r/2 and r] = + 7r/2, for directions 
perpendicular to u and v respectively, we obtain: 



1 

V2 



cos{6 ■ 




(26) 




\^{e,r])) = -^sm{(l>-e) \ COS0 |®( sinc6 |, (27) 



|*(7,(A)) = ^sin(6l-0) I sine I ® I cos*;!. ). (28) 



V2 








Denoting the photon detection in a particular direction with 1 and the non-detection with the joint probabilities 
would be: 

pll = 1 cos2 {0-(f>); = i sin2 {cf> - 9); Pfl = \ sin^ ((? - 0) ; Pf^ = \ cos^ (0 ~ <^)- 

pi^i = icos2(0 -,/)'); Pl^^l^xi^^W Pll = \sxii^{0-(\>'Y Pt^ = \co^{0-4>')\ 

Pll = I C0s2 {0' - 4>) ; P\l = i sin^ (0 - 0' ) ; ^21 = \ Sin^ {& - 0) ; P^l = \ cos^ {0' - 0) ; 
Pll = \ cos^e' - ct>'); Pl^ = \ sin2(</>' - &)■ P% = \ sm\0' - <^')- P°^ = \ cos\0' - <^'). 

If we represent — </>, (f>' — and (j) — 0' hy x, y and z respectively then 0' — (p' , denoted by w, would be —(x + y + z). 

For X = 247.46°, y = 67.49° and z = 157.50° the value of the second term on the left hand side of eq. (23), predicted 
by quantum theory, is | — l/2sin^a; + l/2cos^y + 1/2 sin^^; + l/2cos^ w| = 0.207 and if quantum theory is local, 
according to eq. (23), for the first term we must have 

\Pli - P12 + P2I + P22 I = 0.793. (29) 

However, the value of the loft hand side of the above equality predicted by quantum theory is 1 + 1/2 sin^ x — 
1/2 cos^ y + 1/2 cos^ z + 1/2 sin^ w\ = 1.207. Clearly equality (29) is violated by quantum theory by a factor of 1.522 
which exceeds that of CH results by 0.108. 
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In the presence of white noise the density matrix is 

P = IPnoise + (1 - 7)PQM (30) 

where pnoise = (1/4) x 1 and pqm = |^)(^|- Here 7 is the Bell expression's tolerance of white noise i.e. the maximum 
fraction of white noise admixture for which a Bell expression stops being violated. The joint probability for the above 
state is 

= J + (1 - 7)^6, (31) 

and it is easily scon that the tolerance of equality (29) is 0.586 fraction of white noise which is twice that of CH 
inequality. This is due to the fact that the tolerance is very sensitive to upper bound which is 1 in CH inequality and 
0.793 in our equality. 

Also using the method mentioned in section III we have obtained three inequalities which correspond to different 
values of c and d in eq. (20) which can be tested directly too. One of them, for c = 1 and d = 0, is Clauser-Horne 
inequality, which reads 

< +PA° + P?2 + P21 - P22 < 1- (32) 

The other, for c = 1 and = 1, is 

-1 < +2PA° + ^{'2° - Pl^ - P^? + P^l - P'2^ - P^i < 1- (33) 
Finally, for c = 2 and d = 1, we found 

_1 < _pOO + poi + pio + poo _ pio + pii ^ pOi _ 2pn _ pOO + ^pw < ^ (34) 

A direct proof of this inequality is shown in Appendix B. In terms of x, y, z and w, the inequalities (32), (33) and 
(34) would become 

< 2 (^™^ ^ y ~ w) <1, (35) 

— 1 < sin^ ~\ ~ '"^) ^ 1) (36) 

1 3 

-1 < - + 2 (sin^ ^ ~ sin^ V ~ ^os^ ^ + ^in^ ^ 2, (37) 

respectively. 

The value of inequalities (35), (36) and (37) for x = -67.50°(-37r/8), y = 202.50°(97r/8) and z = -67.50°(-37r/8) 
are +1.207, +1.414 and +2.621 respectively; but for x = 22.50°(7r/8), y = -67.50°(-37r/8) and z = 22.50°(7r/8) 
these are —0.207, —1.414 and —1.621 respectively. While the difference between upper bound and lower bound for 
Clauser-Horne inequality is 1.414, in agreement with the result obtained in [6], in eq. (33) the lower bound itself is 
violated by a factor of 1.414 and the lower bound of eq. (34) is violated by a factor of 1.621 which arc greater than 
that obtained by others for 22(g)22 case (see e.g. [5] and [6]). Please note that the difference between the upper bound 
and the lower bound in the two latter cases are 2.818 and 4.242 respectively and that the range of violation for the 
upper bound for inequalities (32), (33) and (34) are 0.207, 0.414 and 0.621 respectively. These arc the maximum 
values that we could reach and despite the greater range of violation, inequalities (32), (33) and (34) tolerate 0.293 
of white noise which is the same as CH inequality. 

IV. CONCLUSION 

In this paper we have used a numerical method to derive all possible Bell type inequalities in a two party experiment 
namely mm' ®nn' . The dimension of joint probabilities space obtained by this method is not in agreement with the 
results obtained previously in the literature. However we confirmed our numerical results using an analytical method. 
This was done by adopting the assumption of CH, i.e. eq. (8). And on this basis we showed that normalization of 
joint probabilities in a run of experiment and no-signaling in local theories are implied by eq. (1) which is statistically 
trivial. It is an open question that if any assumption other than eq. (8) could leads to these correct results. 
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Wc showed that besides inequahties there are cquahtics that can be used to test local theories. In 22® 22 case 
the equality derived here is violated by a factor of 1.522 tolerating 0.586 fraction of white noise admixture which is 
twice that of the previous restults in the literature. However as in higher dimensions there are inequalities which are 
violated by a stronger factor (see [8]), we expect there may exist equalities in higher dimensions as well which are 
violated by a stronger factor. 

Among many possible Bell type inequalities that we derived, we could only find three types in the form of eqs. (32), 
(33) and (34) which are violated with the same settings used in [5] and [6]. However there are some other settings 
too for which these inequalities are violated and for all of them the inequality in the form of eq.( 34) is violated with 
a stronger factor of 1.621 tolerating 0.293 fraction of white noise admixture. 
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APPENDIX A: DIRECT PROOF OF EQUALITY (23) 



For 22(g) 22 case, using the definition of joint probability P in terms of q's, i.e. 
which are: 



eq. (2), there are 16 possible P's 
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Let's write eq. (23) as 



where Ei and E2 are 



E = \E\\ + LB2 



(Al) 



and 



El - Pll - P12 + P21 + P22 



301 



Eo 



(A2) 
(A3) 



Replacing P's in terms of g's we get 



-£'1 = +9l212 + 9l212 + 9l212 + 9l212 
~9l212 ~ 9l212 " 9l212 ~ 9l212 
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However as g's are all non-ncgativc numbers, Ei 

E = \Ei\ + \E2\ 



\Ei\ and E2 = \E2\. So for local theories we may write 



+9l212 + 91212 + 9l212 + 9l212 + 9l212 + 9l212 + 9l212 + 9l212 



+9l212 + 91212 



9l212 
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91212 



91212 + 91212 



where the last equality is implied by eq. (1) and this justifies eq. (23). 
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APPENDIX B: DIRECT PROOF OF INEQUALITY (34) 



With a similar procedure used in Appendix A, for inequality (34) we have: 
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, r,/ 0100 , 0110 , 1100 , lllON 
+^l9l212 1 9l212 + 9l212 I 9l2127 

— ~(9l212 + 9l212 + 9l212 + 9l212 + 9l212 + 9l212 + 9l212 + 9l212) 

+2(9l212 + 9l212 + 9l212 + 9l212 + 9l212 + 9l212 + 9l212 + 9l212) 

Since according to eq. (1) each of the parentheses on the right hand-side of the above equation is less than 1 and 
greater than zero, it is easily seen that — 1 < J < 2 which is inequality (34). 
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